Post Main Sequence Evolution 


Since a star’s luminosity on the main sequence does not change 
much, we can estimate its main-sequence lifetime from simple 
timescale arguments, and the mass-luminosity relation. If £ œ 
M”, then 


M = 
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If we adopt a main-sequence lifetime of the Sun of 10t? years, 
then ; 
M n 
Ts = 10'° | —— ears 
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Since 7 ~ 3.5, the main-sequence lifetime of a star is a strong 
function of its mass. 


When the mass fraction of hydrogen in a stellar core declines to 
X ~ 0.05 (point 2 on the evolutionary track), the main-sequence 
phase has ended, and the star begins to undergo rapid changes. 


e First, the entire star begins to contract. The core energy gener- 
ation for this stage remains approximately constant, but the star 
increases its luminosity due to the conversion of gravitational en- 
ergy to thermal energy. Simultaneously, the smaller stellar radius 
translates into a hotter effective temperature. For higher mass 
stars, the mass fraction of the convective core begins to shrink 
rapidly. On the diagram, the star travels from point 2 to 3. At 
point 3, the mass fraction of hydrogen in the core is ~ 1%. 


e When hydrogen fusion in a convective core stops, the core’s left- 
over temperature gradient causes energy in the center to flow out- 
ward. Consequently, the core may temporarily cool, and begin 
to contract rapidly, turning some gravitational energy to ther- 
mal energy. The combination of these two effects makes the core 


isothermal. This condition slows the core’s contraction, since the 
demand for energy flow from the center is decreased. (In stars 
with M > 10Mo, the core temperature will never actually de- 
crease, while in stars with radiative cores, core contraction won’t 
actually occur, since it has been adjusting to the decrease in 
hydrogen all along.) 


e While this is happening, the hydrogen rich material around 
the core is drawn inward and eventually ignites in a thick shell, 
containing ~ 5% of the star’s mass. Much of the energy from 
shell burning then goes into pushing matter away in both direc- 
tions. As a result, the luminosity of the star does not increase; 
instead the outer part of the star expands. This “thick shell” 
phase continues with the shell moving outward in mass, until 
the core contains ~ 10% of the stellar mass (point 4). This is 
the Schénberg-Chandrasekhar limit. Stars with larger (by mass 
fraction) cores will reach this point faster than stars with small 
cores. ) 


e Light elements such as lithium, boron, and beryllium fuse at 
temperatures much lower than those of the central core. As a 
result, by the time a star has moved off the main sequence, these 
elements have been destroyed over the inner 98% of the star. 
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The evolution of stars in the HR diagram. For low-mass stars, 
the tracks only extend to the helium flash. 
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The Hertzsprung Gap and the Subgiant Phase 


e When the Schonberg-Chandrasekhar limit is reached, the star 
must change its structure. First, core contraction begins to occur 
on the Kelvin-Helmholtz timescale, and the rapid increase in core 
density causes an increase in the temperatures and densities in 
the shell surrounding the core. The result is a strong increase 
in the rate of nuclear reactions in the shell, which again pushes 
matter away on both sides. The result is a very thin zone of 
hydrogen burning; as the star evolves from point 4 to point 5, 
the mass of the shell will decrease from ~ 3% to ~ 0.5% the mass 
of the star. 


e Although the reaction rates in the shell increase as the shell 
narrows (due to the strong temperature dependence of nuclear 
reaction rates), this is compensated for by the decrease in the 
shell mass. As a result, the total luminosity produced by the 
shell decreases slightly. Moreover, much of this energy goes into 
mechanical work; because the specific heat of the star is nega- 
tive, the increased temperature causes the envelope to expand 
and cool. About half the energy needed to expand the envelope 
comes from shell burning; the other half comes from the envelope 
itself, as it adjusts to the new conditions. The envelope cooling 
rapidly moves the star to the right in the HR diagram, across the 
“Hertzsprung Gap,” until it reaches the limiting Hayashi line. 


e The effective temperature where the star begins to ascend the 
giant branch is approximately independent of mass. It is, how- 
ever, very sensitive to metallicity, since the electrons provided 
by metals are providing most of the surface opacity via H~ ab- 
sorption. Note also that because the convective envelope is large, 
an uncertainty in the mixing length translates into a significant 


displacement of position of the giant branch, with 


ô log Tag ~ |0.02 log = + 0.143 | ó log a ~ 0.16 log a 
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e While in the shell burning phase, the rapidly increasing central 
densities and temperatures can temporarily provide an additional 
source of nuclear energy for the star. While the core of the star 
is mostly helium, there are some trace amounts of metals, in 
particular CNO. Because CNO burning has gone to equilibrium, 
most of the CNO nuclei will be in the form of 4N. As the 
core contracts and the conditions becomes more extreme, 14N 
can fuse with helium, via the reaction 4N(4He, y)!8F(87, v)!8O. 
This rapidly changes all the '*N in the star to 48O, and gives 
the star a little extra life. The place where this occurs depends 
on the mass of the star. Stars with M > 9Mgo fuse 14N in the 
Hertzsprung Gap; as a result, they increase their luminosity on 
their trip across the HR diagram. Stars with M ~ 5Mo burn 
14N on the giant branch. For stars with M < 3Mo, 14N burning 
on the giant branch actually causes them to temporarily descend 
the branch, as the core-shell region adjusts to the new energy. 
None of these effects is directly observable, but 14N will increase 
the lifetime of a subgiant star slightly. 


The Giant Branch 


On the giant branch, the size of the hydrogen-burning shell contin- 
ues to decrease. However, unlike the thick-burning shell phase, the 
decrease in mass is more than offset by the higher nuclear reaction 
rates. To appreciate this, consider that the luminosity produced by 
shell burning depends almost exclusively on the mass of the core. 
The reason for this comes from the equation of hydrostatic equi- 
librium: at the surface of a compact, highly condensed core 


dP GM 


dM Anr4 


Thus the pressure drops by many orders of magnitude in just a 
small region. The extended envelope above the shell is virtually 
weightless, and has almost no influence on the properties of the 
shell. 


The dependence of luminosity on core mass can be estimated from 
homology. If we assume the density, temperature, pressure, and 
luminosity of a burning shell go as some power of the mass and 
radius of the core, i.e., 
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(21.2.1) 


and use homology on the ideal gas law (which is a good approxi- 
mation for regions outside the core), then 


Euta a (Re) 
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giving 
T =+% and h=¢2+% (2122) 


Similarly, the Eulerian equation of hydrostatic equilibrium yields 
dP dr dP'\ GM: Re dP'\ | GM: ; 
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which gives 
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When integrated to a region outside the shell (with negligible pres- 
sure), this equation and (21.2.2) imply 


Tı = Qı +] and T2 = P2 —] (21273) 


In Eulerian coordinates, the luminosity equation is 


dL dr dL! Re dL! 
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If we substitute using € = e9p*T”, then 
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Again, we can integrate over the shell until energy generation van- 
ishes, and re-write the density and temperature in terms of the 
core mass and radius to get 


01 = (A +1)¢1 + wry and 09 = (A+ 1)d2 +yov+3 (21.2.4) 


Finally, we can write the equation for the radiative energy trans- 
port that will take place in the region immediately adjacent to the 
shell. 
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Substituting for k using K = ko P°T", we get 


POOOR 


This can once again be integrated to a region with negligible tem- 
perature, and re-written in terms of the core parameters. The 
result is 


(4 — t)i = 01 + 1 + STI and (4 — t)y2 = o2 + 2 — 1 + ST2 
(21.2.5) 


Equations (21.2.2) - (21.2.5) form a set of 8 (non-linear) equations 
with 8 unknowns. With (a great deal of) effort, these can be solved 
to give 
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Note how the luminosity of a shell burning star behaves. For elec- 
tron scattering (s = t = 0) and CNO burning (A = 1, v = 13), 
La MER. 16/3. for electron scattering and helium burning Ces 


v = 40), L x M} Ro 23/2 Unlike core burning, the luminosity of 
a shell burning star depends both on the core mass and the en- 
ergy generation mechanism. Note also that the temperature of the 
shell has a very simple relation reflective of the virial theorum: 


T x Me/Re. 


We can make further progress by relating the mass of the stellar 
core to its radius. If the core were fully degenerate, this would 
be easy: we could use the mass-radius relation for white dwarfs. 
However, the region immediately below the shell is not degenerate 
(due to the high temperature of the shell), and, although this re- 
gion can have a negligibly small mass, it can occupy a substantial 
fraction of the core’s volume. Thus, a better approximation for the 
core mass-core radius relation of a shell burning star is a two zone 
model consisting of a degenerate core (obeying a polytropic mass- 
radius relation), and an isothermal, ideal gas transition region. 


For the transition region, we can combine the ideal gas law with 
hydrostatic equilibrium to get 


kT dp GM 


um,dr r? 
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This can be integrated from the top of the core, Re, to the place 
where electron degeneracy takes over, r;. The result is 


In p — ln pe = - (= ) [2 — i} (2122.7) 


This equation can quickly be simplified. First, since the transition 
region is at the same temperature as the shell, then from (21.2.6), 


T x Me/Re, and the term in parenthesis is a constant (C1). Next, 
we can estimate p; by equating the ideal gas and electron degen- 
eracy density laws. From (7.3.8), we have 


3 3 
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Similarly, we can substitute for rą using the polytropic mass-radius 
relation for white dwarfs 


n = a [EET eRT M 


= 9.71 x 102 u75 M53 = Cupp 5M7" (16.1.4) 
and for ln pe using the homology relations for (21.2.6) 
In pe = C3 + Q1 ln Me + doln Re 


This gives us a relation between Re, T, u, He, and Me 
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(21.2.9) 
For reasonable values, this equation yields 
dln Re 
~ —0.1 212 4 
dln Me ne ( 0) 


(a rather weak dependence). The net result is that 


L x MŽ (21.2.11) 


with z ~ 8 for CNO burning shells, and z ~ 15 for helium burning 
shells. 


Once we have an equation for the luminosity of a shell burning star 
as a function of core-mass, it is then relatively easy to calculate 
the luminosity of the star as a function of time. Consider that as 
a Shell burning star evolves, it continually deposits more mass on 
its central core, with 
dM. £ 
dt XQ 


where X is the mass fraction of the fuel, and Q the amount of 
energy generated by one gram of material. Thus 


ti 
f da M-7dM (21.2.12) 
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where t = 0 is the fiducial time when £ = Lo, and K is the constant 
of proportionality relating £L to Me. Equation (21.2.12) is easily 
integrated to yield 


QX >% = 
aan Me, aM Jat 


which, after a bit of manipulation, becomes 
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This is a sharply increasing exponential! The star’s evolution ac- 
celerates dramatically as it ascends the giant branch. 
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